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Abstract 

The paper is concerned with a problem of coherent (measurement-free) filtering for physically realizable 
(PR) linear quantum plants. The state variables of such systems satisfy canonical commutation relations and are 
governed by linear quantum stochastic differential equations, dynamically equivalent to those of an open quantum 
harmonic oscillator. The problem is to design another PR quantum system, connected unilaterally to the output 
of the plant and playing the role of a quantum filter, so as to minimize a mean square discrepancy between the 
dynamic variables of the plant and the output of the filter. This coherent quantum filtering (CQF) formulation is 
a simplified feedback-free version of the coherent quantum LQG control problem which remains open despite 
recent studies. The CQF problem is transformed into a constrained covariance control problem which is treated 
by using the Frechet differentiation of an appropriate Lagrange function with respect to the matrices of the filter. 

£ ; I. Introduction 

Interconnection of open systems, whose internal dynamics are affected by interaction with the 
surroundings, is often engineered so as to stabilize the resulting network of such systems via 
redistribution and dissipation of energy generated by active nodes or coming from the environment. In 
addition to its role in redirecting the energy flow, interaction provides a universal mechanism for creating 
correlations whereby current states of different subsystems acquire and store dynamic "footprints" of 
^! each other and of the past history of the whole system. This informational aspect of interaction is directly 
employed in the classical Kalman filter whose state is continuously updated by the measurement process 
from a stochastic system, thus enabling the filter to develop and maintain relatively strong correlation 
qh. with the unknown state of the system. The ability of such a filter to track a classical linear system 
1—1 with finitely many degrees of freedom and known dynamics by extracting as much information from 
the noisy observations as possible is, in principle, limited only by digital implementation. The situation 
is qualitatively different in regard to estimating the dynamic variables of a quantum stochastic system 
which are noncommutative operators on a Hilbert space evolving in time according to the laws of 
quantum mechanics [11]. The measurements, which result from the interaction of the quantum system 
with a relatively invasive classical device, are accompanied by irreversible loss of quantum information 
as a consequence of the projection postulate of quantum mechanics [5]. It is the idea of using a 
system of the same kind, that is, another quantum system, weakly coupled to the quantum mechanical 
object of interest (a quantum plant), which underlies the coherent (that is, measurement-free) quantum 
j> ■ filtering/control paradigm. This approach replaces measurement with interaction of quantum systems, 
possibly mediated by light fields, where the energy flow can be employed for stabilization/control [7] and 
the quantum information manifests itself through quantum correlations between the dynamic variables 
of the systems in the course of time. An important class of quantum stochastic systems is provided 
by open quantum harmonic oscillators [2], [3] whose variables satisfy canonical commutation relations 
(CCRs) and are governed by linear quantum stochastic differential equations (QSDEs) driven by boson 
fields [14]. In combination with the preservation of CCRs under unitary evolutions in the Heisenberg 
picture of quantum dynamics, the specific energetics of such systems (quadratic Hamiltonian and linear 
coupling to the external fields) imposes certain constraints [8], [16] on the coefficients of the governing 
linear QSDEs in order for them to be physically realizable (PR) as open quantum harmonic oscillators. 
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Such systems can be implemented in practice by using quantum-optical components [3]. Being quadratic 
with respect to the state-space matrices, the PR constraints make the coherent quantum counterpart [13] 
of the classical LQG control problem [9] substantially harder to solve than its classical predecessor. 
In fact, the coherent quantum LQG (CQLQG) feedback design problem remains open despite recent 
studies [18], [19] which explore different approaches to its solution. In the present paper, we consider an 
infinite-horizon coherent quantum filtering (CQF) problem for PR linear quantum plants. The question 
of interest is to design another PR quantum system, connected unilaterally to the output of the plant and 
playing the role of a quantum filter, so as to minimize a steady-state mean square discrepancy between 
the dynamic variables of the plant and the output of the filter. In the absence of measurements and in 
the presence of PR constraints, the machinery of recursive Bayesian estimation (including conditional 
expectations), so useful for Kalman filtering, is inapplicable to the CQF problem. Following [18] based 
on algebraic ideas from [1], [17], we transform the CQF problem into a constrained covariance control 
problem which is treated by using the Frechet differentiation of an appropriate Lagrange function with 
respect to the matrices of the filter. Since the CQF setting is a simplified feedback-free version of the 
CQLQG control problem mentioned above, this leads to a more explicit set of algebraic equations for 
the state-space matrices of an optimal PR quantum filter which are amenable to further analysis to 
be published elsewhere. Note that the recognition of the need to take into account PR constraints in 
coherent quantum filtering problems dates back to [4], although that work was mainly concerned with 
measurement-based mean square optimal filtering. A coherent quantum filtering problem has recently 
been discussed in [12], where, unlike the present paper, optimization of the filter was not considered. 

II. Physically Realizable Quantum Plant 

As in the linear quantum control settings [8], [13], [18] mentioned above, the quantum plant considered 
below is an open quantum stochastic system with canonically commuting variables whose internal 
dynamics are affected by the quantum noise from the environment. More precisely, the plant has n 
dynamic variables x\ (t),.. -,x n (t) which are self-adjoint operators on the tensor product Hilbert spaceQ 
J£f (g> J^i evolving in time t and satisfying the canonical commutation relations (CCRs) 

[ XiX T ] ■= ([xj,x k ]) 1<jik<n = xx T - (xx T ) T = 2i&\. (1) 

Here, x := (xk)i^k^n is the vectoid of the plant variables (the time argument is often omitted for the 
sake of brevity), [r\,Q := rjC, — C,r\ is the commutator of operators, i := \J— 1 is the imaginary unit, 
and 0i is a constant real antisymmetric matrix of order n (the subspace of such matrices is denoted 
by A„) which is assumed to be nonsingular, so that n is even. The operators xi(0), ...,x n (0) act on the 
initial complex separable Hilbert space Jt\ of the system, and &\ is the boson Fock space [14], which 
provides a domain for the action of the quantum Wiener processes w\ (?),..., w mi (t). The latter are 
self-adjoint operators on &\ (which are obtained from pairs of field annihilation and creation operators 
by using complex unitary (2 x 2) -matrices) and represent the quantum noise from the environment with 
the quantum Ito table dwdw T = f^df. Here, w := (wjt)i<jfc<mi» and is a constant complex positive 
semi-definite Hermitian matrix of order m\. Similarly to (OQ), the imaginary part J\ := Im^i G A mi of 
the quantum Ito matrix Q.\ specifies the CCRs between the quantum Wiener processes as 

[dw,dw T ] =2iJidt. (2) 

In what follows, the real part Re^i is the identity matrix I m of order ni\, so that 

Qi=/ mi +*7i. (3) 

'Some of the spaces and parameters associated with the plant are equipped with the subscript "1", whereas the subscript "2" is used 
for analogous objects pertaining to the quantum filter introduced in Section Hill 

2 Vectors are organized as columns unless specified otherwise, and the transpose (-) T acts on matrices with operator-valued entries as 
if the latter were scalars. 
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Also, it is assumed that the noise dimension m\ is even, and the CCR matrix 7i has a canonical form 



l 

o 



Mi 



'Mi 



'Ml 





(4) 



where <g> is the Kronecker product of matrices, and \i\ :=m\/2. The plant state vector x evolves in time 
and contributes to a ^-dimensional output of the plant y (whose entries are also self-adjoint operators 
on Jt[ <g)#i) according to QSDEs 



dx = Axdt + Bdw, dy = Cxdt + Ddw. 



(5) 



Here, A e R nxn , B e R nxmi , C G W xn , D e W xmi are given constant matrices, with A Hurwitz. In 
addition to the asymptotic stability, the plant is assumed to be physically realizable (PR) as an open 
quantum harmonic oscillator [2], [3] whose Hamiltonian is quadratic and the coupling operators are 
linear with respect to the plant variables. By the results of [8], [13], [16], in the case of linear quantum 
dynamics being considered, the PR property is equivalent to the algebraic relations 

A0i + 0iA T + fi/ifi T 



0iC T +5JiD T 



0, 




(6) 
(7) 



which describe the preservation in time of the CCRs between the state and output variables of the plant 
described by (OQ) and [x,y T ] = 0. Indeed, ©, © are obtained from the relationships 

d[x,x T ] = (A[x,x T ] + [x,x T ]A T + 2iBJiB T )dt, 



d[x,y T ] = (A[x,y T ] + [x,x J }C J + 2iBJ 1 D T )dt 



(8) 
(9) 



which follow from the bilinearity of the commutator [11], the quantum Ito rule, and the commutativity 
between adapted processes and forward increments of the quantum Wiener process [14], in view of ©, 
©. 

III. Physically Realizable Coherent Quantum Filter 

Consider another PR open quantum stochastic system with q canonically commuting dynamic 
variables which are influenced by the environment and the plant in a unilateral fashion. Due to this 
cascade connection, which is shown in Fig. [T] the state and output variables of the second system under 




Fig. 1. The dynamics of the filter are influenced by the quantum Wiener process CO and the plant output y through a unilateral cascade 
connection described by (|5j, Jilt . 



consideration acquire quantum correlation with the plant variables in the course of time, which enables 
this system to be regarded as a coherent (that is, measurement-free) quantum filter. A performance 
criterion for such a filter will be specified in Section ITVl Now, the filter state variables £i (f ),..., 
are assumed to be self-adjoint operators on the Hilbert space J#\ <g> J^2 <8> &\ <8> &i satisfying the CCRs 

[<^ T ]=2z0 2 , (10) 

where £ := (E,k)i^k^q an d ©2 £ A^, with q even and det©2 7^ 0. Here, M2 is the initial space of the 
filter and is the boson Fock space which provides a domain for the action of a quantum Wiener 
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process ft) := (ft)/t)iscfcscm 2 - The latter commutes with the plant noise w and drives the filter variables 
according to the QSDE 



dt, = a£,dt + bdco + edy, drj = ct,dt + ddco 1 



(11) 



where a G R qxq , b G R qxm2 , c G W xq , d G W xm \ e G R qxp are constant matrices, and a is Hurwitz. 
The quantum Wiener process ft) is assumed to have even dimension ni2 and a canonical quantum Ito 
table 



dcodco = Q. 2 dt, 



^2 := Imi+iJi 



analogous to ©, ©. Here, J% G A OT2 is a CCR matrix of the filter noise in the sense that 

[dft>,dft) T ] = 2iJ 2 dt, h :-- 



I, 
-hi ' 



(12) 



(13) 



with /i2 := mj/l. The matrices Z?, e in (fTTT) play the role of gain matrices of the quantum filter with 
respect to the filter noise ft) and the plant output y. The plant, filter and the environment form a closed 
quantum system which is governed by ©, (fTTT ). Therefore, the (n + q) -dimensional vector 



St :-- 



formed by the plant and filter state variables, is driven by the combined quantum Wiener process 



(14) 

(15) 

(16) 
(17) 



The matrices sd , £3 have a block lower triangular structure due to the absence of feedback, as opposed 
to the closed-loop quantum control settings [8], [13], [18]. The plant and filter noises w, ft) result from 
interaction of the systems with external boson fields which are assumed to be in the product vacuum 
state v := Vi <E> t>2 on the space &\ (g) J^2- Since the noises commute with each other and are uncorrected, 
then, in view of ©-© and (fl~2l) . (fT3~l) . the combined Wiener process W in (fl~5T) has a block diagonal 
quantum Ito table 



of dimension m := m\ +111,2 according to the QSDE 

dS£ = £rfS£dt + mW, 

with 



A 


38\= 


B 





eC a 




eD 


b 



dWdW v =Qdf, 



Hi 
i2 2 



Here, J G A m denotes the corresponding CCR matrix of 

[dW,dW T ] =2iJdt, 



Ji 
h 



(18) 



(19) 



In a similar vein, the plant and filter variables are assumed to commute, so that the combined vector 
(fl~4l) has a block-diagonal CCR matrix: 



\SC, eT T ] = 2/0, 0:= 



@i 

02 



(20) 



4 



where use is made of (Q~K (flOT ). Due to the unitary evolution of the isolated system formed by the plant, 
filter and their environment, the CCR matrix is also preserved in time, which is equivalent to the 
algebraic Lyapunov equation 

£/® + ®£/ T + £gJ^ T = 0. (21) 

The left-hand side of (12TI) is a real antisymmetric matrix whose diagonal block, the upper off-diagonal 
block and the second diagonal block are computed as 

A0i + 0iA T + BJ X B T = 0, (22) 
(&iC T + BJ l D T )e T = 0, (23) 
a@ 2 + ®ia T + eDJiD T e T + bJ 2 b T = (24) 

in view of (fTTT) . (fT9l) . (1201) . The fulfillment of d22j), ([23]) is guaranteed by the PR properties ©, © of 
the quantum plant for an arbitrary coherent quantum filter (1TTI) . whereas (|24l) and the equality 

® 2 c T + bJ 2 d T = (25) 

describe PR conditions for the filter which correspond to the preservation of the CCRs (flOT ) and [£, , r] T ] = 
0, respectively. The derivation of the PR conditions (|24l . (|25l) is similar to that of ©, © in ([8]), ®, 
except that the dynamic variables of the filter are driven not only by the quantum noise ft), but also by 
the plant noise w through the plant output y according to the QSDE (PT61) . which explains the presence of 
the additional term eDJ\D T e T in (|24l) . Since 02 is nonsingular, the general solution of (1241 . considered 
as a linear equation with respect to the matrix a, is described by 

a = 2& 2 R ~ ^(eDJiD T e T + bJ 2 b T )® 2 1 . (26) 

Here, R is an arbitrary real symmetric matrix of order q (the subspace of such matrices is denoted by 
E) q ) which specifies the quadratic Hamiltonian <^ T i?<^/2 of an equivalent representation of the PR filter 
as an open quantum harmonic oscillator. Furthermore, (1251 allows the matrix c to be expressed in terms 
of b as 

c = -dJ 2 b T & 2 1 . (27) 

The coupling of the state-output matrix c to the noise gain matrix b makes the optimization of the 
coherent quantum filter (flTT ) qualitatively different from that of the classical controllers and filters, 
regardless of the performance criterion. Indeed, (l27l shows that the PR quantum filter requires an 
"intake" of the additional quantum noise ft) (through b ^ 0) in order to produce a useful output r\ with 
a nonzero "signal" component 

(28) 

In their original form [8], [13], [16], the PR conditions also involve a specification of the noise 
feedthrough matrices D, d in ©, (fTTT) as those formed from rows of orthogonal matrices, whereby 

DD T = I p , dd T =I r . (29) 

Therefore, p ^ m\, r ^ m 2 , and both D and d have full row rank. If y in ([5]) were a classical observation 
process, the full row rank property of D would correspond to nondegeneracy of the measurements. Also, 
since d is of full row rank, the map M. qxm2 3 b^ c E W xq , described by (l27l . is surjective, so that the 
matrix c can be assigned any value by an appropriate choice of b. Although (|29l will simplify algebraic 
manipulations, it is the rank properties of D, d that are principal for what follows. 
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IV. Coherent Quantum Filtering Problem 

Consider the problem of constructing a PR coherent quantum filter (flTT) of fixed dimensions (and 
with a fixed noise feedthrough matrix d), described in Section |TTT1 so as to minimize a steady-state mean 
square discrepancy between the filter output variables and the state variables of a given PR quantum 
plant © specified in Section HO More precisely, let Sf{t) := denote an s-dimensional 

quantum process defined by 

2?:=Fx-GC, =^JT, (30) 

where F G t sx ", G G f sxr are given matrices, with s ^ r and G having full column rank (the role of 
this assumption is clarified later), and £ is the signal part (12~8l) of the filter output r\ in (flTT) . so that, in 
view of (fill) . 

«*f : = [F -Gc] . (31) 

The coherent quantum filtering (CQF) problem is formulated as the minimization of the quantity 

S := lim E(3f 7 3?) ll = i^ T ^,P) ll 

— ^min, subject to dill),©. (32) 

Here, the 1/2 factor is introduced for further convenience, Sf 1 = Yfk=i ls me sum °f squared 
entries of iF, and E(-) denotes the quantum expectation over the product state (D®V, where 07 is the 
initial quantum state of the plant-filter composite system on J4?\ (g) J^2, and V is the vacuum state of the 
external fields. Also, (M,N) :=Tr(M*N) is the Frobenius inner product of complex or real matrices, 
with (•)* := ((-)) T the complex conjugate transpose, and P is the real part of the steady-state quantum 
covariance matrix of the vector S£ in (fl4l) : 

K:= lim E(^T^T T ) =P + i®, P:=ReK, (33) 

with lim^ +00 E^ = 0, since the matrix srf in (fTTT) is Hurwitz. The latter condition is ensured by the 
Hurwitz property of the matrix a, since stf is block lower triangular and A is Hurwitz. The matrix K 
has the imaginary part Im&T = given by (1201) if the coherent quantum filter is also PR. In view of 
(fT6l) . (fi~8T) . the matrix P in (1331) is a unique solution of the algebraic Lyapunov equation 

£/P + P£/ T + 3§3§ T = (34) 

and coincides with the controllability Gramian [9] of the pair £§). The minimum in the CQF 
problem d32) is taken over the quadruple (a,^,c,e) G R qxq x Ri xm 2 x W xq x R qx P of the state-space 
matrices of the filter (flTT ) subject to the PR constraints (1241 . (1231) . with a fixed noise feedthrough matrix 
d E W xm2 satisfying (l29l . In particular, if F = G = I n , with r = s = n, the CQF problem (1321 consists 
in approximating the plant state vector x by the signal part £ of the filter output r\ from (1281 ) so as to 
minimize the "estimation error" x — C, in the mean square sense. 

V. Qualitative Dependence on Filter Matrices 

The performance index in the CQF problem (|32l is a composite function (a, c, e) n- , ^) t?, 
where the triple G M( n +«) x («+«) x R( n +i) xm x M sx ("+^) of matrices from (HI), (|3B depends 

affinely on the matrix quadruple (a,Z?,c,e) (with both being regarded as elements of appropriate direct 
sum Hilbert spaces). The smoothness of <§ with respect to sd ', £$, c € (regardless of the specific structure 
of these matrices) is ensured by the smooth dependence of the controllability Gramian 

P = A^(S33f) (35) 

from (l34l) on srf ', £$ over the open subset of Hurwitz matrices srf. Here, A a denotes the inverse Lyapunov 
operator, that is, a particular case of the inverse Sylvester operator E a p , which is associated with Hurwitz 
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matrices a, /3 and maps an appropriately dimensioned matrix X to the unique solution Y :=L a p(X) 
of the algebraic Sylvester equation aY + Yj3 +X = 0: 

A a :=E aaT , Z a)P (X) := J^o^X^'dt. (36) 



For what follows, the controllability Gramian P in (|35l) (and other related matrices of order n + q) 
is split into blocks P\\ e §„, P22 € § 9 , Pi 2 = Pji e K" X9 and the corresponding block-rows Py, and 
block-columns P,£ in accordance with their association with the plant and filter variables in (fl4l as 



Pu Pn 

P21 P22 



Pu 
P2. 



[P.i P.2]. (37) 



The Lyapunov equation (1341) . whose left-hand side is a symmetric matrix with a similar block 
partitioning, can be written in terms of (l37T i as 

A/ 5 n+/ 5 nA T + J B5 T = 0, (38) 

AP l2 + Pi2a T + PnC T e T +BD T e T = 0, (39) 

aP 22 + Pz2<* T + eCPi2 + P 2 iC T e T + ee T + bb T = 0, (40) 

where (fTTT ). (|29l are used. The inverse Lyapunov and Sylvester operators (1361 allow the solution of 
(|38T>-(|40l) to be represented as 

P n =A A (BB T ), (41) 

P n = Z AM T((P n C T + BD T )e T ) 1 (42) 

P 22 = K{eCP n + Pj 2 C T e T + ee T + bb T ). (43) 

Since the matrix Pn in (I4TT) is specified completely by the quantum plant © and does not depend on 
the coherent quantum filter, the matrix P\\C +BD T in (l42l) is also constant. Therefore, for any given 
Hurwitz matrix a of the filter (fTTT) . the matrix P12 is a linear function of e. Hence, P22 in (1431) is a 
homogeneous quadratic polynomial of the filter matrices b, e from (ITTb . whose coefficients depend on 
a and which does not contain the bj^e^ cross-terms with mixed entries of b, e. The performance index 
§ in the CQF problem (1321 is representable in terms of the block partitioning (|37l) as 

S= (F T F,P n )/2-(F T Gc,P l 2) + (c 7 G T Gc : P 2 2)/2, (44) 

where (I3TT) is used together with the symmetry of the matrices and P. In combination with the 

quadratic dependence of P22 on b, e, and the linear dependence of P12 on e discussed above, (1441) 
implies that $ is quadratic in b, c, e for any fixed Hurwitz matrix a. Such dependence of the quadratic 
performance index S on the filter matrices also holds in the classical filtering problem. However, in 
contrast to its classical predecessor, the CQF problem (1321 is constrained by the PR conditions (1241) . 
(1251) which couple the skew-Hamiltonian part of the matrix a in (1261) to b, e and the matrix c in (1271) 
to fc. 

VI. Conditions of Stationarity 

If an optimal filter exists for the CQF problem (|32l . such a filter is among stationary points of the 
Lagrange function 

if : =S + (S, a& 2 + & 2 a T + eDJ x D T e T + bJ 2 b T ) /2 

+ (r, 2 c T +^ T ) 

=<f - (S0 2 , a) + (S, eDJiD T e T + bJ 2 b T )/2 

-(rj7 2 ,^) + (r T 2 ,c). (45) 
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Here, the 1/2 factor is introduced for further convenience, and S G A q , T e M, qxr are Lagrange multipliers 
associated with the PR constraints (|24l) (whose left-hand side is a real antisymmetric matrix) and (I25T ). 
respectively. We will be concerned with a quadruple (a,Z?,c,e) of unconstrained state-space matrices 
of the coherent quantum filter (fTTT) which is a stationary point of the Lagrange function jSf in (|45l) for 
given matrices S, T, with a Hurwitz. The Lagrange multipliers S, T are to be found so as to ensure that 
the filter (which depends parametrically on S, T) satisfies the PR conditions. In order to find stationary 
points of the Lagrange function Jzf in (143T) . we will compute its Frechet derivatives by using the chain 
rule and the following lemma [18, Lemma 2] based on algebraic techniques from [1], [17]. 

Lemma 1: The Frechet derivatives of the function <§ in (l32l) with respect to the matrices si ', ^ 
from (fTTT ). (l3"TT ). with srf Hurwitz, are computed as 

d jg £ = QP=:H, d®g = Qm, d^S = ^P, (46) 

where P is the controllability Gramian from (|34l . and Q := A ia /r( < jf 1 V) is the observability Gramian 
of the pair (s/^) satisfying the algebraic Lyapunov equation 

^ T Q + Q^ + ^ T( tf = 0. (47) 

□ 

The matrix H in (1461) will be referred to as the Hankelian of the triple (stf : 33, ^) since the spectrum 
of H is formed by the squared Hankel singular values [9], [17] of an appropriate linear time-invariant 
system. A block-wise form of (1471 is given by 

A T <2i i + QnA + C T e T Q 2l + Qi 2 eC + F T F = 0, (48) 
a T Q2i+Q2iA + Q 2 2eC-c T G T F = 0, (49) 
a J Q 2 2 + Qua + c T G T Gc = 0, (50) 

which is similar to (I38l)-(l40l) except that the lower off-diagonal block is considered instead of the upper 
one. The solution of (|48l - (|50l) is found by using the inverse Lyapunov and Sylvester operators (I36T ) as 

fin =A A T(C T e T Q 2 i + Qj l eC + F T F) : (51) 
<22i = ^ a T A {Q 22 eC - c T G T F), (52) 
fi 22 = A fl x (c T G T Gc), (53) 



which corresponds to (|4Tb — (143D . The block fi22 in (|53l is computed first and is then substituted into 
(1521) in order to find Q 2 i, while fin in (I5TT) is irrelevant for further discussions. 

Lemma 2: The Frechet derivatives of the Lagrange function .if from (1431) with respect to the state- 
space matrices a, b, c, e of the quantum filter (flTT) . with a Hurwitz, are computed as 

d a J? = H 22 -Z® 2 , (54) 

d b ^ = Q 22 b - EbJ 2 - TdJ 2 , (55) 

d c ££ = -G T FP l2 + G T GcP 22 + r T 2 , (56) 

d e ^ = H 2l C T + Q 2 iBD J + Q 22 e - ZeDJ { D J : (57) 

where use is made of the block partitioning of the controllability and observability Gramians P, Q and 
the Hankelian H according to (T3~7l . □ 
The proof of Lemma |2] is given in Appendix |A] The following conditions of stationarity of the 
Lagrange function Jzf are obtained by equating its Frechet derivatives to zero. 
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Theorem 1: The coherent quantum filter (fTTT) . with a Hurwitz, is a stationary point of the CQF 
problem (1321 if and only if there exist Lagrange multipliers S G A q , T e IR ?xr in (1451) such that the 
equalities 

#22 = S0 2 , (58) 

Q 22 b = ZbJ 2 +TdJ 2 , (59) 

G T GcP 22 - G T FP l2 = -r T 2 , (60) 

^2iC T + <22i5D T + <2 2 2e = SeD7 1 D T (61) 

are satisfied together with the PR conditions (1241 ). (|23T > . □ 
Together with the Lyapunov equations (1341) . (1471) and the PR conditions (1241) . (1251) . the relations (|58T>- 
(I6TI) form a complete set of algebraic equations for finding an optimal filter among stationary points 
(a,b,c,e,E,T) of the Lagrange function in (|45l) . Since (l59~l) . (l60l) can, in principle, be solved for b, c, 
the Lagrange multiplier T can be found so as to satisfy the PR condition (1251) . Also, (16TT) can be solved 
for e, as discussed below. In view of (1581) . which describes the stationarity of the Lagrange function 
with respect to a, and the antisymmetry of S, the matrix H 22 = 02 1 2 S0 2 is skew-Hamiltonian in the 
sense of the symplectic structure specified by ® 2 Note that (1581) is implicit in a. In the next section, 
we will obtain a more explicit equation and discuss its solvability with respect to a along with that of 
(EUhED) for b, c, e. 

VII. Combining the Stationarity and Lyapunov Equations 

We will need two identities for the Gramians P, Q from (1341) . (1471) and the Hankelian H from (|46T ) 
which hold regardless of whether the QSDEs ©, (fTTT ) are PR. To this end, we introduce a matrix 

T := Q 2 .£?P. 2 = Q 2 iAP 12 + Q 22 eCP 12 + Q 22 aP 22 , (62) 

where use is made of (flTT ) along with the block partitioning of P, Q, H in accordance with (1371) . 
Lemma 3: The matrix T, defined by (1621) . with ^ Hurwitz, satisfies 

T + (H 2l C T + Q 2l BD J + Q 22 e)e T 

+H 22 a T + Q 22 bb T = 0, (63) 

T + a T H 22 + c T G T (GcP 22 - FP n ) = 0. (64) 

□ 

The proof of Lemma |3] is given in Appendix |B] We will now combine the general identities (1631) . 
(1641) with the equations (158l)-(l6"TT) of stationarity of the Lagrange function Jzf from Theorem Q] and the 
PR conditions (|24l . (l25l) . 

Lemma 4: If the PR coherent quantum filter is a stationary point of the Lagrange function in (|43T) . 
with a Hurwitz, then the matrix Ea + Tc is symmetric, and the matrix T from (|62l satisfies 

T=(Sa + rc)0 2 . (65) 

□ 

The proof of Lemma |4] is provided in Appendix O The lemma will be used to obtain a more 
explicit equation for a than (|58l) . Now, suppose the "filter blocks" of the controllability and observability 
Gramians P, Q are both nonsingular: 

^22 y 0, £22 y 0. (66) 

The positive semi-definiteness P 22 y 0, <2 2 2 y is inherited from P, Q. By the general form of 
Heisenberg's uncertainty principle [5], the steady-state quantum covariance matrix of the filter state 
variables satisfies 

P 22 + i& 2 = lim E(^ T ) ^ 0, (67) 
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L 


= Q 2 iQ2u 


M 


-P12P22, 


N 


— Q22 


S 


= ®2P 2 2i 


T 


= e£r, 


u 


= (22202^22 



which is a stronger property than P22 ^ 0. The assumptions (1661) enable the following matrices to be 
defined: 

(68) 
(69) 
(70) 

Here, only N, T involve the Lagrange multipliers S, T from (I45T ). while L, M, S, U are completely 
specified by 02 and the appropriate blocks of the Gramians P, Q, with neither P\\ nor Q\\ being 
involved. The eigenvalues of the matrix S in (|69l ) are purely imaginary and symmetric about the origin, 
with the spectral radius satisfying 

r(5)<l (71) 

in view of (l67l) . The matrix "ratios" from (168T>— (TTOl) allow the equations of Theorem [TJ to be made more 
explicit in a, b, c, e. 

Lemma 5: Under the assumptions (1661) . the stationarity equations (l58l) -(l6"TT) for the CQF problem 
(1321) can be represented in terms of the matrices L, M, N, S, T , U from (|68T>— (|70T> as 

LM = NS-I q , (72) 
b = NbJ 2 + TdJ 2l (73) 
c=(G T Gy\G T FM-T T U), (74) 
e = NeA-(L(P n C T + BD T )+P 2 iC T ) 1 (75) 

where 

A:=D7iD T . (76) 
Furthermore, the matrix a of such a filter satisfies the relation 

a = (Na + Tc)S-(LA + eC)M. (77) 

□ 

The proof of Lemma [5] is given in Appendix ID] Note that (1741) employs the full column rank assumption 
on the matrix G in (|30l) to ensure the invertibility of G T G. While (1741) is already solved for c, we will 
show how the stationarity equations (fTTb . (1731) . (1751) can be solved for the other filter matrices a, b, e. 
To this end, let [[[oc,/3]]] : X \-t aXfi denote the linear operator of the left and right multiplication of 
an appropriately dimensioned matrix X by given real matrices a, /3, respectively. More generally [18, 
Section 7], for any positive integer g and compatibly dimensioned matrices a\,...,a g and jSi, . . . , j8„, a 
special linear operator of grade g is defined by 

[[[ai,j8i|...|a g , J 8J:=£[[[a,, J 8J, (78) 

where the matrix pairs are separated by "|"s. If, in each of the pairs, the matrices a^, are either both 
symmetric or both antisymmetric, then the operator ([78} is self-adjoint with respect to the Frobenius 
inner product of matrices. Now, the solvability of (1771 ) with respect to a depends on whether the grade 
two special linear operator 

IV, I -N, 5]]] = IQ22 1 , P22I [[[222, ^22 I -S, 02]]] (79) 
(with o denoting the composition) is invertible on Mfl xq . In the case of invertibility, ([771) takes the form 

« = Wq, h I ~N, 5]]] -\TcS-(LA + eC)M) . (80) 
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In a similar vein, (1731) . (1751) can be solved for b, e, respectively, as 

b=lI q ,I m \-NM]~\TdJ 2 ), (81) 
e = -gl q ,l p | -N,AT l (L(P n C T + BD T ) + /> 21 C T ), (82) 

provided the following special linear operators of grade two are invertible: 

IkM I -NM = flfc&V-J o iQnM I -S, J 2 ffl, (83) 
p 9 ,7 p | -N, A]]] = [[[e^ 1 , / J o [[[ 222 , / p | -S, A]]] . (84) 

By combining the spectral property (1711 of the matrix S from (1691 ) with similar properties 

r(7 2 ) = l, r(A)<l (85) 

of the matrices 7 2 m (fl3l) and A in (|76l) . and applying Lemma [6] of Appendix |E] to the operators in (|79|) . 
(l83l . (184b . it follows that the condition 

r(N) < 1 (86) 

is sufficient for the invertibility of all three operators. Note that, similarly to [18, Proof of Lemma 
5], the second of the spectral relations ([85) follows from © and (T2~9l) which imply that the complex 
Hermitian matrix zA = iDJ{D T satisfies —I p =4 *A ^ I p . Under the assumption Q 22 y 0, the condition (l86l) 
is equivalent to the strict convexity of the Lagrange function Jz? in (1451) with respect to b. More precisely, 
Jzf inherits quadratic dependence on b from S, as discussed in Section [V] and, in view of (I55T) . the 
second order Frechet derivative c^=£? = [[[<222,7 m2 | — S,7 2 ]]] is a grade two special self-adjoint operator 
on the Hilbert space Wi xm2 whose positive definiteness is indeed equivalent to (l8~6l) by Lemma |6[ Here, 
we have again used the spectral property of the matrix 7 2 from (I83T ). 

VIII. An Iterative Algorithm Outline 

A reasoning, similar to that in the previous section, shows that, under the assumptions (I66T ) and (l8~6l) . 
the matrices b, c in (1731) . (|74[) satisfy the PR constraint (1251) if and only if the matrix T in (l70l) is related 
to £ by 

r = [/- T (M T F T G + 2 1 Z77 2 J T G T G), (87) 

where (-) _T := ((-) _1 ) T is the composition of the matrix inverse and transpose. This allows T to be 
eliminated from ([811 as 

b=l[I q J m2 \-N,J 2 

| -U T & 2 \j 2 d T G T GdJ 2 ^\U- T M T F T GdJ 2 ) 1 (88) 

which involves an invertible special operator of grade three. The matrix b from ([88) can be substituted 
into (ITTT) in order to find c. In an iterative algorithm for numerical computation of an optimal quantum 
filter in the CQF problem (l32l . the relations (188) . (|8~7l . (|27l can be employed to update the matrices b, 
T, c for given matrices M, N, U from (|68T) - (|70l) . In a similar vein, (f80T >. (l8~2l can be used for updating 
the matrices a, e for given L, M, N, S, T and P\ 2 . This is accompanied by updating the matrices L, M, 
S, U according to (T68T> — (T70T> in terms of the blocks of the Gramians P, Q which are computed for given 
filter matrices a, b, c, e as described by (l42l) . (1431) and (l52l) . (153T ). The algorithm loop can be closed by 
computing the block H 22 of the Hankelian H in (l46l) and updating the Lagrange multiplier S in (l69l) 
as E := (H 22 ®2 1 + ® 2 1 7^ 22 )/2. In the case n = q, the algorithm can be initialized with the state-space 
matrices of a classical Kalman filter from the next section. 
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IX. Reduction to the Classical Kalman Filter 

If the PR constraints (El), ([25]) are made inactive by letting E = 0, T = in <@3]) (so that the CCRs 
become irrelevant), the Lagrange function corresponds to a classical filtering problem. In this case, the 
matrices N, T in (|69l) , (1701) vanish, and, under a simplifying assumption that the filter has the same 
state dimension q = n as the plant, the necessary conditions of optimality (I72l-(l75l). (1771) take the form 



LM = -I n , (89) 

a = -(LA + eC)M, (90) 

b=0, (91) 

c=(G T Gy l G T FM 1 (92) 

e = - L(P n C T + BD T ) - P 2 iC T = -L(nC T + 5D T ), (93) 

where, in view of (1681) . the real positive semi-definite symmetric matrix 

U:=P n -MP 2 i (94) 



is the Schur complement [6] of P22 in (1371) . In particular, (|9Tb shows that the additional noise co 
(uncorrected with the plant noise w) is redundant in the optimal filter, in conformance with the classical 
Kalman filtering theory. Since (l89l) implies that L= — M , an appropriate similarity transformation 
^4ff^,fl4 oao~ l , b 1 — y ob, c h-> c<7 , e 1— >■ ae of the filter (flTT) with a nonsingular matrix o G M nXM 
(which does not have to be symplectic in the sense that o&20 T = ©2) leads to L = —I n , M = I n . The 
corresponding matrices a, c, e in (|90l ), (|92l , (|93l reduce to 

a=A-eC, c = (G T G) l G T F, e = TlC T + BD T , (95) 

with c being specified completely by the matrices F, G. Accordingly, the Lyapunov equation (i38T>— (|40T> 
takes the form of an algebraic Riccati equation [10] for the matrix IT from (l94l) : 

An + ita t + bb t - (nc T + bd t ) (cn + db t ) = o. (96) 

In view of (1291) . both (1931) and (1961) indeed correspond to the classical Kalman filter SDE &£, =A^dt + 
e(dy — C^dt), with dy — C^dt the innovation. 
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Appendix 



A. Proof of Lemma [2] 
Since the matrices 



^ do not depend on a, then, in view of (fTTT) and Lemma [U the first 



variation of the function $ with respect to the matrix a is 8 a <% = ( d^<§ , 



o 

5a 



= (#22) So). Hence, 

d a S = H22, which, in combination with d a (Z&2,a) = Z&2 from (145T ). implies (1541) . Similarly, since 
sd ', ^ do not depend on the matrix b, then, in view of (fTTT) and Lemma [H the first variation 



of the function $ with respect to b takes the form S^S* = (d.ggS, Sb j = (Q2ib,8b), so that 

— Qnb- In combination with d/,(Z 1 bJ2b T ) = —2ZbJ2 and dt,(TdJ2,b) = TdJ2 from (|45T ). this 
establishes (1551) . In a similar vein, since the matrices s^, SS do not depend on c, then, in view of 
(I3TI) and Lemma CD the first variation of the function <§ with respect to the matrix c is computed as 
8 c <g = (dtfS', [0 -G8c}) = - (G Tc &P % 2, 8c) . Therefore, d c S = -G T tfP. 2 = -G T FP i2 + G T GcP 22 , which, 
in combination with <9 c (r T ®2,c) = T T 02 from (l45l ). yields (l56l) . Finally, since ^ in (|3TT) is independent 



of the matrix e, then, in view of (fTTT) and 



to e is 8 e S = ( B^S 





SeC 



+ ( d.^ 



-emma fU the first variation of the function <S with respect 
(H 2 iC T + (Q 2 iB + Q 22 eD)D T ,8e), and hence, 





SeD 



d e S = H 2 iC T + Q 2 iBD r + Q 22 e, where use is made of (|29]>. In view of d e (Z, eDJiD T e T ) = -IZeDJiD 1 



in (1451) . this leads to (1571) . thus completing the proof of the lemma. 



B. Proof of Lemma \3\ 

Left multiplication of both sides of the Lyapunov equation (l34l) by Q2, and right multiplication of 
the Lyapunov equation (|47T) by P,2 yields 



Q 2 .^P + H 2 .^ T + Q 2 .t 



0, 



/T H. 2 + Q^P.i + tf T tfP.2 = 0. 



(Bl) 
(B2) 



Here, use is made of the identities Qj,P = Hj, and QP»k — H.k, which follow from the definition of 
the Hankelian H in (1461) and its block partitioning as in (|37T >. The definition of the matrix in (fTTT) 
implies that 



#2. & = [Hi 1 A J H 2 i C T e J + H 22 a r } 



l H. 2 



A T Hi2 + C T e T H 22 
a T H 22 



(B3) 
(B4) 



By recalling the matrices SB, from §VT§, d2B, and substituting {QS$Sf) 2 2 = {Qi\B + Q 2 2eD)D 1 e T + 
Q 2 2bb T and (^ T ^P) 2 2 = c T G T (GcP 22 - F P n ) together with (|29]>, (|B3]), (Ell), into the appropriate blocks 
of (iBll) . (IB21) . the resulting equations lead to f)3]), (J64J). ■ 
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C. Proof of Lemma |4] 

Substitution of the stationary conditions (l58i (|59i (I6TT) and the PR conditions ([24]), ([25]) into the 
identity (|63l) yields 

T + SeD/iDV + S0 2 a T + ZbJ 2 b T - TcQ 2 
= T + Z(eDJiD T e T + & 2 a T + bJ 2 b T ) - Tc& 2 

= Y-(Sa + rc)02 = O, (CI) 

thus proving (I65T ). In a similar vein, by substituting (1581) . (1601) into (1641) . it follows that 

T + a T Z® 2 - c T r T 2 = T + (a T Z - c T r T )0 2 = 0. (C2) 

Since the left-hand sides of (IC II) . (|C2I) are equal to each other, and det©2 7^ 0, then Za + Tc = c T T T — 
a T Z. The latter equality, in view of E T = — Z, is equivalent to the matrix Za + Tc being symmetric. ■ 

D. Proof of Lemma \5\ 

The equation (1721) is obtained by the left multiplication of both sides of (1581) by Q 22 and right 
multiplication by P 22 , followed by using the identity H 22 = Q 2 \P\ 2 + Q 22 P 22 together with the matrices 
L, M from (|68) and N, S from (J6^]». The equality dT3]) is obtained by the left multiplication of d59]) by 
Q22 and using the matrices N, T from (1691 ), (P70i The equation (1741 is obtained by the left multiplication 
of (1601) by (G T G) _1 (secured by the full column rank of G) and right multiplication by P 22 and using the 
identity T T & 2 P 22 l = (Q 22 l T) T Q 22 & 2 P 22 l = T J U which follows from (|70]>. The equation d75]) is obtained 
by the left multiplication of (16TI) by QI, and using the identity H 2 i = Q 2 iPn + Q 22 P 2 \ along with the 
matrices L, N from (1681) . (169T ). Finally, (|75T) is obtained by the left multiplication of (l65l) by an d 
right multiplication by P 22 and using the representation T = Q 22 (LAM + eCM + a)P 22 for the matrix T 
in d62]) in terms of L, M from (|68]>- ■ 

£. Positive Semi-definiteness of Grade Two Operators 

Lemma 6: Suppose a G S r , j8 G § s and a G A r , T G A s , with a >- 0, /3 >- 0. Then a criterion of positive 
semi-definiteness of the grade two special self-adjoint operator 4> := [[[oc, jS | c,t]]], defined by (1781) on 
the Hilbert space W xs , is given by 

cj>^0 <^=> r(«- 1 c7)r(T/3" 1 ) < 1. (El) 

□ 

Proof: Consider the decomposition of 4> into the sum 4> = ^ + of grade one self-adjoint 
operators := [[[cu,/3]]] and U := [[[c,r]]], where y due to the assumptions on the matrices oc, 
/3. Hence, 4> = V¥(I + A) V^, where A := ^-1/2^-1/2 = ^-1/2^-1/2^-1/2^-1/2]]] is a grade 

one self-adjoint operator whose spectrum is symmetric about the origin [18, Lemma 1]. In view of this 
symmetry, the condition r(A) ^ 1 is not only sufficient, but is also necessary for the fulfillment of 4> y- 0. 
It now remains to note [18, Section 7] that r(A) = r(a _1 / 2 aa _1 / 2 )r(/3 _1 / 2 T/3 -1 / 2 ) = v{a~ l a)r{x^~ l ), 
and the equivalence (lETI) is proved. ■ 
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